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Abstract 



Based on the analysis of a large sample of RXTE/PCA observations of several 
black hole binaries in the low spectral state we show that a correlation exists 
between the spectral parameters and characteristic noise frequency. In particular, 
the amplitude of reflection increases and the slope of Comptonized radiation 
steepens as the noise frequency increases. We consider possible implications of 
these results on the accretion flow models and discuss a possible observational 
test aimed to discriminate between different geometries of the accretion flow. 
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Comptonization of soft seed photons in a hot, optically thin electron cloud 
near the compact object is thought to produce the hard X-ray radiation in the 
low spectral state of accreting black holes (Sunyaev & Truemper 1979, Sunyaev 
& Titarchuk 1980). The slope of the Comptonized spectrum is governed by 
the ratio of the energy deposited into the electrons and the influx of the soft 
radiation into the Comptonization region; the lower the ratio the steeper the 
Comptonized spectrum (e.g. Sunyaev & Titarchuk 1989, Haardt & Maraschi 
1993, Gilfanov et al. 1995). 

Reflection of the Comptonized radiation from neutral or partially ionized 
matter, presumably the optically thick accretion disk, leads to appearance of 
characteristic features in the spectra of X-ray binaries (Fig. 1.1). The main 
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Figure 1.1 The counts spectra of Cyg X-l at different epoch plotted as a ratio to a power law 
model. The power law photon index is indicated against each spectrum. 



signatures of the emission reflected from cold neutral medium are well known 
- the fluorescent K Q line of iron at 6.4 keV, iron K-edge at 7. 1 keV and a broad 
hump above ~ 20 — 30 keV (Basko, Sunyaev & Titarchuk 1974, George & 
Fabian 1991). The exact shape of these spectral features in the X-ray binaries 
depends on the ionization state of the reflecting medium (Ross, Fabian & Young 
1999) and might be modified by the strong gravity effects and intrinsic motions 
in the reflector (e.g. Fabian et al., 1989). The amplitude of the reflection 
signatures depends primarily on the ionization state and on the solid angle of 
the reflector as seen from the source of the primary radiation. 

Based on the analysis of a large sample of Seyfert AGNs and several X-ray 
binaries Zdziarski, Lubinski & Smith 1999 (see also Zdziarski 2000) found a 
correlation between the amplitude of reflection and the slope of the underlying 
power law. They concluded that the existence of such a correlation implies a 
dominant role of the reflecting medium as the source of seed photons to the 
Comptonization region. 

The power density spectra of the black hole binaries in the low state (see van 
der Klis 1995 for a review) plotted in the units of frequency x power usually 
appear as a superposition of two or more rather broad humps containing most of 
the power of aperiodic variations below several tens Hz (Fig. 1.2). These humps 
define several characteristic noise frequencies of which the lowest usually 
corresponds to the break frequency v\, T of the band limited noise component. 
Although a number of theoretical models was proposed to explain the power 
spectra of X-ray binaries (e.g. Alpar et al. 1992, Stella Vietri & Morsink 
1999), the nature of the characteristic noise frequencies is still unclear. Despite 
the fact that the characteristic noise frequencies vary from source to source and 
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Figure 1.2 Left: The power density spectra of Cyg X-l at different epoch. The power density 
spectra are plotted as frequency x (power density), i.e. in units of Hz x rms 2 /Hz. The 
observations and the symbols are the same as in Fig. 1.1 Right: The relation between the QPO 
centroid frequency and the break frequency of the band limited noise for Cyg X-l, GX339^1 
and GS 1354-644. The solid line is f qpo = 6 x f br . 



from epoch to epoch, they are correlated with the break frequency of the band 
limited noise (Fig. 1.2, Wijnands & van der Klis 1999). 

Gilfanov, Churazov & Revnivtsev 1999 showed that a tight correlation exists 
between the characteristic noise frequencies and the spectral parameters. Based 
on RXTE/PCA observations of Cyg X-l from 1996-1998 they found that 
the amplitude of reflection increases and the spectrum of primary radiation 
steepens as the noise frequency increases. An identical behavior was found for 
another black hole candidate, GX339-4, by Revnivtsev, Gilfanov & Churazov 
2000. Interpreting these correlations Gilfanov, Churazov & Revnivtsev 1999 
suggested that increase of the noise frequency might be caused by the shift of 
the inner boundary of the optically thick accretion disk towards the compact 
object. The related increase of the solid angle, subtended by the disk, and of 
the influx of the soft photons to the Comptonization region lead to an increase 
of the amount of reflection and steepening of the Comptonized spectrum. 

Below we investigate further the relation between the spectral parameters 
and characteristic noise frequencies based on a larger sample of objects and 
observations and discuss implications of these results on the geometry of the 
accretion flow in the low spectral state of black hole binaries. 
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1. THE DATA 

Observations and data reduction. The results presented below are based 
on the publicly available data of Cyg X-l, GX339^1 and GS 1354-644 obser- 
vations with the Proportional Counter Array aboard the Rossi X-ray Timing 
Explorer (Brandt Rotschild & Swank 1996) from 1996-1998 during the low 
(hard) spectral state of the sources. In total our sample included ss 60 obser- 
vations of these sources. The energy and power density spectra were averaged 
for each individual observation. The energy spectra were extracted from the 
"Standard Mode 2" data and ARF and RMF were constructed using standard 
RXTE FTOOLS v.4.2 tasks. The "VLE" and "Q6" models were used for the 
background calculation with the preference being given to the "VLE" model 
when possible. A uniform systematic error of 0.5% was added quadratically to 
the statistical error in each energy channel. The power spectra were calculated 
in the 2-16 keV energy band and in the « 0.005 — 32 Hz frequency range fol- 
lowing the standard X-ray timing technique and normalized to units of squared 
fractional rms. 

The spectral model. The energy spectra were fit with a model consisting of a 
power law with a superimposed reflected continuum (pexrav model in XSPEC, 
Magdziarz & Zdziarski 1995) and an intrinsically narrow line at ~ 6.4 — 6.7 
keV. The centroid energy of the line was a free parameter of the fit. For 
all sources the binary system inclination was fixed at i = 50 degrees, the 
iron abundance - at the solar value of Ap c = 3.3 • 10~ 5 and the low energy 
absorption - at Nu = 6 • 10 21 cm~ 2 . The effects of ionization were not 
included. In order to approximately include in the model the smearing of the 
reflection features due to motion in the accretion flow the reflection continuum 
and line were convolved with a Gaussian, which width was a free parameter 
of the fit. The spectra were fitted in the 4-20 keV energy range. This spectral 
model is identical to that applied by Gilfanov, Churazov & Revnivtsev 1999 to 
the Cyg X-l data, except that in the present analysis we let the line centroid 
energy be a free parameter of the fit. Strength of the reflected component is 
characterized in the pexrav model by the reflection scaling factor R, which 
approximately measures the solid angle subtended by the reflecting media as 
seen from the source of the primary radiation, R ss Q r efi / 2ir, so that R = 1 
for an isotropic point source above an infinite optically thick slab. 

The spectral model is obviously oversimplified. Moreover, several important 
parameters, such as the binary system inclination angle and the iron abundance 
were fixed at fiducial values. Therefore the best fit values do not necessarily 
represent the exact values of the physically interesting parameters. Particularly 
subject to the uncertainties due to the choice of the spectral model is the 
reflection scaling factor R ~ Q/2ir. This might explain the values of R 
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Figure 1.3 The best fit values of the photon index of the underlying power law plotted vs. 
reflection scaling factor for Cyg X-l, GX339-4 and GS 1354—644. The solid line shows the 
dependence T(R) of the photon index of the Comptonized radiation T on the strength of the 
reflected emission R expected in the disk-spheroid model assuming the disk albedo a = 0.1, 
Thomson optical depth of the cloud tt = 1 and the ratio of the temperature of the seed photons 
to the electron temperature Ttb /T e = 10 -4 . The dashed line shows T(R) dependence expected 
in the plasma ejection model (Beloborodov 1999) for a = 0.15, r T = 1, T bb /T e = 3 ■ 10~ 3 
and (is = 0.3. 



exceeding unity obtained for some of the spectra. However, as was shown by 
Gilfanov, Churazov & Revnivtsev 1999 the model correctly ranks the spectra 
according to the strength of the reflected component and the slope of the 
underlying power law. 

Power spectra approximation. The power density spectra in units of rms 2 / Hz 
were fit in the 0.005-15 Hz frequency range with a model, consisting of two 
broken constant components and a Lorentzian profile P v = i + ( 2 (u-u )//iu)' 1 • 
We use centroid frequency of the Lorentzian profile as a characteristic noise 
frequency. Its dependence on the break frequency of the band limited noise 
(the first broken constant component in the model) is shown in Fig. 1 .2. 
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Figure 1.4 The QPO centroid frequency plotted vs. reflection scaling factor for the same set 
of observations of Cyg X-l, GX339-4 and GS 1354-644 as in Fig. 1.3 and 1.6. 



Correlation between the spectral parameters and the noise frequency. 

The results of the energy and power spectra approximation are presented in 
Fig. 1.3 and 1.4 showing dependence of the slope of the primary emission and 
characteristic noise frequency on the amplitude of the reflected component. 
The noise frequency and the spectral parameters are clearly correlated. An 
increase of the amplitude of the reflected component is accompanied with a 
steepening of the spectrum of the primary emission and an increase of the noise 
frequency 

2. DISCUSSION 

Slope of the Comptonized spectrum and reflection. The photon index T 
of the Comptonized radiation is governed by the Compton amplification factor 
A which equals to the ratio of the energy deposited into the hot electrons to 
the influx of the soft seed photons to the Comptonization region. The concrete 
shape of the T(A) relation depends on the ratio /T e of the temperatures of the 
seed photons and the electrons, the Thomson optical depth and the geometry. In 
the simplest although not unique scenario the correlation between the amplitude 
of the reflected component R and the slope of the primary radiation T (Fig. 1.3) 
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Figure 1.5 Relation between the Compton amplification factor and photon index of Comp- 
tonized radiation for an isotropic point source of the soft seed photons with temperature T r in 
the center of a uniform sphere with electron temperature T e and Thomson optical depth tt = 1 
(hollow symbols) and tt = 2 (solid symbols). The symbols are a result of the Monte-Carlo 
calculations, the solid curves were calculated according to the formula given in the text. 

could be understood assuming that there is a positive correlation between the 
solid angle subtended by the reflecting media, $l re fi, and the influx of the soft 
photons to the Comptonization region (Zdziarski, Lubinski & Smith 1999). 
Existence of such correlation is a strong argument in favor the reflecting media 
being the primary source of the soft seed photons to the Comptonization region. 
In the absence of strong beaming effects a correlation between Q re fi and the 
seed photons flux should be expected since an increase of the solid angle of the 
disk seen by the hot electrons (= O re /0 should generally lead to the increase 
of the fraction of the disk emission reaching the Comptonization region. 

In order to demonstrate that such geometrical effects can explain the observed 
dependence F(R) we consider two idealized models having different cause of 
change of the O re //. In the first, disk-spheroid model, an optically thin uniform 
hot sphere with radius R sp h is surrounded with an optically thick cold disk with 
an inner radius Rdisk (e.g. Poutanen Krolik & Ryde 1997), the Q r efi depending 
on the ratio Rdisk / Rsph- Propagation of the disk towards/inwards the hot sphere 
(decrease of Rdisk/ Rsph) leads to an increase of the reflection scaling factor, 
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a decrease of the the Compton amplification factor A and a steepening of 
the Comptonized spectrum. In such context the model was first studied by 
Zdziarski, Lubinski & Smith 1999 and we used their results to calculate the 
relation between the reflection scaling factor R and the Compton amplification 
A. In the second, plasma ejection model, value of the tt re fi is defined by the 
intrinsic properties of the emitting hot plasma, particularly by its bulk motion 
with mildly relativistic velocity towards or away from the disk, which itself 
remains unchanged (Beloborodov 1999). In the case of an infinite disk, values 
of reflection below and above unity correspond to the hot plasma moving 
respectively away and towards the disk. Both models predict relation between 
reflection R and Compton amplification factor A which can be translated to 
T(R) given a dependence T(A) of the photon index of Comptonized spectrum 
on the amplification factor. The latter was approximated by: 



This formula is based on a representation of the Comptonized spectrum by a 
power law in the energy range 3fcT^ — 3kT e and agrees with the results of the 
Monte-Carlo calculations with reasonable accuracy for optical depth tt ~ 1 
and T hb /T e ~ 1CT 5 - 10" 3 (Fig. 1.5). 

The expected T{R) relations are shown in Fig. 1.3. With a proper tuning 
of the parameters both models can reproduce the observed shape of the T(R) 
dependence and in this respect are virtually indistinguishable. The observed 
range of the reflection R ~ . 3 — 1 and the slope F ~ 1.5 — 2 . 2 can be explained 
assuming variation of the disk radius in from Rdisk ~ Rsph to Rdisk ~ in 
the disk-spheroid model or variation of the bulk motion velocity from v ~ 0.4c 
away from the disk to v ~ in the plasma ejection model. Finally, it should be 
emphasized that these idealized models do not include a number of important 
effects which might affect the particular shape of the T(R) dependence. 

Relativistic smearing and reflection. The spectrum of emission, reflected 
from a Keplerian accretion disk is expected to be modified by special and 
general relativity effects (e.g. Fabian et al., 1989). If increase of the reflection 
is caused by a shift of the inner boundary of the accretion disk closer to the 
compact object, a correlation between the reflection and relativistic smearing 
should be expected. Such a correlation is a generic prediction of the model and 
might be used to discriminate between different assumptions about geometry 
of the accretion flow. The energy resolution of RXTE/PCA is not adequate 
to study relativistic smearing of the reflection features with sufficient degree 
of confidence. However, the data shown in Fig. 1.6 might present an evidence 
in favor of correlated behavior of the reflection and relativistic smearing. If 
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Figure 1.6 The best fit value of the relative width of the relativistic smearing of the reflection 
features plotted vs. reflection scaling factor for Cyg X-l, GX339-4 and GS1354-644. 



confirmed, it would prove that the spectral and timing properties in the low 
state of the black hole binaries are defined mainly by the position of the inner 
boundary of the optically thick accretion disk. 

3. CONCLUSIONS 

In the low spectral state of black hole candidates the amplitude of reflection 
R, slope of the Comptonized radiation T and characteristic noise frequency 
change in a correlated way. The correlation between the slope and reflection 
indicates that the accretion disk is the primary source of the soft seed photons 
for Comptonization. We considered two idealized models, the disk-spheroid 
and plasma ejection model, and showed that with appropriate tuning of the 
parameters either model can explain the observed F(R) dependence. The 
RXTE data also hints at a correlated change of the width of relativistic smearing 
of the reflection features. If confirmed by observations with better energy 
resolution, this correlation may help to distinguish between different geometries 
of the accretion flow. 
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